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Abstract 

We discuss alternative iteration methods for differential equations. We provide a 
convergence proof for exactly solvable examples and show more convenient formulas 
for nontrivial problems. 



1 Introduction 



The asymptotic iteration method (AIM) is an algorithm for the exact and 

approximate solution of second-order ordinary differential equations [1, 2]. 
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It has been applied to a wide variety of problems that include exact and 
approximate calculations of eigenvalues and eigenf unctions, implementation 
of perturbation theory, nonrelativistic and relativistic problems, among others 



Almost all those references are devoted to applications of the approach and 
little has been done to provide a sound foundation for the AIM. The proofs 
given in the first two papers on the subject [1,2] do not account for all the 
properties of the method. Recently, Matamala et al derived a most interesting 
connection between the AIM and continued fractions [12]. 

The purpose of this paper is to provide a deeper insight into the AIM. In Sec. 2 
we give well known formal solutions to the differential equation. In Sec. 3 we 
discuss the AIM and some conditions for its successful apphcation to exactly 
solvable and nontrivial problems. We discuss a simple problem with constant 
coefficients in order to illustrate the kind of solutions expected from the AIM. 
In Sec. 4 we propose alternative iterative approaches that we apply to the same 
simple problem just mentioned. In Sec. 5 we apply the AIM to the harmonic 
and anharmonic oscillators, compare alternative recurrence relations and the 
effect of different initial conditions. Finally we summarize our results and draw 
conclusions in Sec. 6. 



2 Differential equation 

The purpose of this paper is the exact or approximate solution of differential 
equations of the form 



[3-11]. 
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where p{x) and q{x) are integrable functions. We can easily obtain a formal 
solution to this equation if we factorize it as 



d , , 
— + a[x) 
ax 



y{x) = 0. 



(2) 



On comparing both equations we realize that a{x) 
is a solution of the Riccati equation 



—p{x) — b{x) where b{x) 



h'{x) — hixY — p{x)h{x) + q{x) = 0. 



(3) 



Straightforward integration of equation (2) yields the general solution 



y{x) = exp 



X t XL 

- j h{u) du lCi + C2 j exp j {2b{u) + p{u)) du 



dt\ (4) 



where Ci and C2 are integration constants. Notice that it is a linear combina- 
tion of two independent solutions of Eq. (1). 



The logarithmic derivative 



y{x) 



(5) 



also satisfies the Riccati equation (3): 



f{x) - f{xf - p{x)f{x) + q{x) ^Q. 



(6) 



All those results are well known, and we summarize them here merely to 
facilitate the discussion below. 
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3 Asymptotic Iteration Method 



The AIM is applicable if p{x) and q{x) are C°° functions. If we differentiate 
Eq. (1) n times we obtain 



(7) 



where 



Pn=Pn-l + PPn-1 + Qn-l, 
qn = qn-l + QPn-l, ^= 1,2, 

Po^p, Qo = q- 



(8) 



Ciftci et al [1] proved that if 



On _ Qn-l 
Pn Pn-1 



a 



(9) 



then 



y{x) = exp 



a{u) du 



X i 

C1 + C2 j exp j {2a{u) +p{u)) du 



dt } .(10) 



is a general solution of the differential equation (1). It is clear that Eq. (10) 
agrees with Eq. (4) if a{x) — b{x) [2] and its existence is therefore independent 
of the condition (9). Moreover, we conclude that a{x) should satisfy the Riccati 
equation (3). 

Saad et al [13] proved that the exact solutions just discussed are polynomial 
functions. In what follows we provide an alternative proof that is more con- 
venient for the treatment of nontrivial problems. If wc differentiate the ratio 
Qn-i/Pn-i a-nd use the recurrence relations (8) we obtain 



Qn-l 



kPu-I. 



Qn-l , 

-P 



Pn-1 



pI-1 
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where 



K = ?nPn-l - qn-lPn- (12) 

Therefore, if 5„ = and Pn-i 7^ then a = qn-i/Pn-i satisfies the Riccati 
equation (3) and y = C exp[— a{u) du] is a solution to the differential 
equation (1). Since — Pn-i{y' + ocy) = we conclude that y{x) is 

polynomial of degree at most n . Conversely, if y{x) is a polynomial solution 
of degree n, then y' + Qn-iy/Pn-i = 0, provided that ^ 0, a = Qn-i/pn-i 
satisfies the Riccati equation and 5„ = 0. Summarizing, there is a polynomial 
solution y{x) to Eq. (1) if and only if 5„ = and Pn-i 7^ 0. This is exactly 
theorem 2 of reference [13] except that present proof does do not require that 
Pn 7^ 0. One can easily verify that if p„ = then Qk — Pk — ^ for all /c > n 
under the conditions above. 

The condition (9), although useful for exactly solvable problems, is not suitable 
for nontrivial ones where we require that 

lim — — a. (13) 

n— >oo r) 
J/n 

Notice that if 

lim 4^ = (14) 

then a{x) given by Eq. (13) is a solution of the Riccati equation (3) and Eq. 
(10) is the general solution of the differential equation (1). 

In order to understand some of the main features of the AIM it is convenient 
to discuss a simple problem already considered earlier. If p{x) and q{x) are 
constant, then the AIM recurrence relations (8) are exactly solvable: 

Qn^qPn-l, 
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Pi,2 = ^, A = Vp2 + 4g, (15) 

where the constants Ci and C2 are determined by the conditions po — p and 
P-i = 1- 

If |pi| < IP2I equation (13) gives us a = — pi that is a root of the Riccati 
equation (3) (b'{x) = 0). We appreciate that the AIM yields the root with 
smaller modulus and we then obtain the general solution by means of Eq. 
(10). This extremely simple exactly solvable example is interesting because its 
solutions are not polynomials. 

If |pi| = IP2I and pi 7^ P2 the AIM does not converge but we can overcome 
this difficulty quite easily. The function v{x) — y{x) exp(/3x) is a solution 
of the differential equation v"{x) = pv\x) + qv{x), where p = 2l3 + p and 
q = q — pf3 — (3'^ . The new roots arc pi,2 = (3 + pi,2 and the AIM converges 
for the modified differential equation because |pi| 7^ |p2|. We will see that 
transformations of this sort are useful for the treatment of the Schrodinger 
equation. Clearly, this strategy fails when pi — p2- 



4 Iterative Riccati method 



It is clear from the discussion above that the Riccati equation (6) is central 
to the AIM. We can derive the AIM from the Riccati equation if we look for 
a solution of the form 

/(.) ^ (10) 

On substituting this expression into the Riccati equation and rearranging con- 
veniently we obtain [2] 

B B' + A + pB ^ ' 
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If we solve the equations A — A' + qB, and B — B' + A + pB iteratively we 
obtain the AIM recurrence relations 



An — A^_i + qBn-i, 

E„ = 5;_i + A,_i+p5,_i, (18) 

except that we do not have any prescription for the initial conditions and 
Bq. However, according to the example of the preceding section, the initial 
conditions do not appear to be that relevant for nonpolynomial solutions. The 
most important fact is that 

/ An-l \' f An-l \^ p ^'^~^ \ q— (19) 

where (5„ = B^-iA^ — B^A^-i, which clearly tell us that the sequence of ratios 
An-i/ Bn-i may converge to a solution of the Riccati equation. 

We can rearrange the equation for A and 5 in a different way 
A A' -pA + qB 



B B' + A 



(20) 



and derive the alternative recurrence relations 



An = A^„i - pAn-l + qBn-1 

Bn = B'^_, + An-i. (21) 

This approach does not agree with the AIM. We appreciate that the iterative 
Riccati method is somewhat more arbitrary than the AIM. 

The sequences {A^p , B^p} and {A^p , B^p} given by the recurrence relations 
(18) and (21), respectively, are not independent. In fact, it is not difficult to 
prove that A^^\x) = A(/)(x)e"(^), and B^^^x) = S^i)(a;)e"(^), where u'{x) = 
p{x). Both recurrence relations should give the same result if the initial con- 
ditions are also related by the same transformation. 
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If we apply equations (18) to the exactly solvable example discussed above we 
obtain the AIM result. On the other hand, if we apply equations (21) the result 
is the root of the Riccati equation with greater modulus: lim„^oo(74„/i?„) = p2- 
In both cases we assume A'q = B'q = 0. 

5 The Schrodinger equation 

Direct application of the AIM to the Schrodinger equation 



may lead to divergent sequences for nonpolynomial solutions. Notice that we 
meet the same difficulty when we apply the AIM to the simple example above 
with p = 0. In order to overcome it we make the transformation ip[x) ~ 
g{x)y{x) and apply the AIM to the resulting differential equation for y{x): 



It has been shown that the rate of convergence of the AIM sequences depends 
on the function g{x) [2]. For example, in the case of the harmonic oscillator 
V{x) — x^ we choose g{x) — exp(— x^/2) and y{x) satisfies the Hermite dif- 
ferential equation that we discuss briefiy below. Notice that in the case of 
eigenvalue problems one has to determine the value of the energy E together 
with the solution y[x). We obtain the eigenvalues E from the roots of 8^ = 0. 
In the case of polynomial solutions this equation yields exact eigenvalues for 
finite n, but for nontrivial problems we obtain increasingly accurate results as 
n^oo. [1,2] 

In order to have a deeper insight into the approaches derived above it is conve- 
nient to consider a simple differential equation with variable coefficients p{x) 



ijj"(x) = [V{x) - E]ilj{x) 



(22) 




(23) 
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and q{x). One of the simplest examples is the Hermite differential equation [14] 

y"{x) = 2xy'{x) - 2my{x) (24) 

with polynomial solutions for m = 0, 1, This example is different from the 

preceding one in that the iteration method determines the value of m and, 
consequently, of q together with the solutions of the corresponding Riccati 
equation. 

The terminating condition yields the values of m corresponding to the Hermite 
polynomials as shown by: 

S„ = -2"+^m(m - 1) . . . (m - n) = 0. 

n = l,2,... (25) 
Besides, each of the functions 



m = 0=^-^ = 0, n>0 

Bn 

An 1 

m = 1 =^ — = — , n > 

Bn X 

An Ax 

m — 2=^ — = , n> 1 

Bn 2a;2-l' - 

„2 



An 3(2a;^ - 1) 

m — — = — n> 2 

Bn x(2x2-3)' - 



(26) 



satisfies the Riccati equation with the corresponding value oiq — —2m. In this 
case we have chosen the AIM initial conditions Aq = q, and Bq = p. Notice 
that given the value of m we obtain the exact polynomial solution from An/ Bn 
for any n > m — 1. 

If, on the other hand, we choose, for example, Aq = Bq = 1, then, for a given 
m we obtain the corresponding Hermite polynomial for all n > m + 1. That 
is to say, we need more iterations for the same result, which suggests that the 
AIM prescription is most convenient for this case. 
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As suggested by the simple example in Sec. 3 the alternative recurrence rela- 
tion (21) may yield other kind of solutions. We have confirmed this point in 
the case of the Hermite equation. The roots of 5„ = are negative integers 
m = —1, —2, . . ., and An/ Bn = —2x, —(1 + 2x^)/x, . . . for sufficiently large 
but finite n are solutions to the corresponding Riccati equation. 

However, if we choose the initial conditions Ao{x) = q{x)e^^ , and Bq = p{x)e^^ 
in the recurrence relations (21), then we obtain exactly the AIM results in 
complete accordance with the discussion in Sec. 4. 

As a nontrivial model we consider the anharmonic oscillator V{x) — x^. We 
follow an earlier application of the AIM and choose g{x) — exp(/3x^/2) where 
P is an adjustable parameter [2]. However, in this case we select the initial 
conditions Aq = 1 and Bq = 1 arbitrarily for the recurrence relations (18). 

Results are similar to those given by the standard AIM [2] which shows that 
the initial conditions are not so relevant in the case of nonpolynomial prob- 
lems. Particularly, if a great number of iterations is required as in the present 
application. 

6 Conclusions 

In this paper we try to provide an alternative proof for the AIM in the case of 
exactly solvable examples and develop equations that appear to be more con- 
venient for the discussion of nontrivial problems. We also expect to place the 
AIM in a more general context of iterative algorithms for differential equa- 
tions as Matamala et al have also done regarding the continued fractions 
algorithm [12]. 

Present results suggest that the AIM gives a convenient prescription for the 
starting point of the recurrence relations if one is looking for the square inte- 
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grable solutions of the Schrodinger equation. However, other initial conditions 
may lead to identical results, particularly in the case of nonpolynomial prob- 
lems where a great number of iterations is necessary for accurate results. 

One can derive alternative recurrence relations that also give solutions to 
the Riccati equation and, consequently, to the linear differential equation. 
However, those alternative recurrence relations require that one chooses the 
initial conditions carefully; otherwise one may obtain unwanted solutions. 
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